Delocalization in Coupled Luttinger Liquids with Impurities 
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We study effects of quencfied disorder on coupled two-dimensional arrays of Luttinger liquids (LL) 
as a model for stripes in high-Tc compounds. In the framework of a renormalization-group analysis, 
we find that weak inter-LL charge-density-wave couplings are always irrelevant as opposed to the 
pure system. By varying either disorder strength, intra- or inter-LL interactions, the system can 
undergo a delocalization transition between an insulator and a novel strongly anisotropic metallic 
state with LL-like transport. This state is characterized by short-ranged charge-density-wave order, 
the superconducting order is quasi long-ranged along the stripes and short-ranged in the transversal 
direction. 

PACS numbers; 71.10.Hf, 71.27.-|-a, 71.23.-k 



I. INTRODUCTION 

Quasi-one-dimensional electron liquids play a paradig- 
matic role in describing the conductive properties of 
variety of physical systems such as organic conductors,! 
quantum-HaJL systemsjj and striped phases in high-T( 
compounds.Ba Recent studies of weakly coupled Lut- 
tinger liquids (LL) have provided evidence for the sta- 
bility of p|ji[|LjFermi-liquid behavior in more than one 
dimensio 
Fermi gas 



as opposed to results for an isotropic 2D 
This remarkable result is a consequence of 
the combined effect of single particle/Cooper pair tunnel- 
ing and Coulomb interactions between the LLs. Earlier 
studies either excluded hopping,El treated aliinter-LL in- 
teractions separately as weak perturbationalj or focused 
on single particle tuaneling for strong repulsive intra- 
LL interactions only.tj In Refs. it was shown that 

backscattering and particle hopping processes between 
the LLs can be irrelevant for sufficiently strong inter- 
LL forward scattering. The resulting state was called 
"sliding Luttinger liquid" (SLL). For a large range of in- 
teractions, these processes can be partially relevant and 
lead to charge-density wave (CDW), trans»erse super- 
conductor (SC) or Fermi Liquid (FL) phases Hu'El Experi- 
ments have paovided evidence for ID transport in high-Tc 
compounds.E3 Theoretically, novel and to date essentially 
unexplored behavior can arise from disorder, which is in- 
duced by doping in these materials. 

Here, we examine the role of electron scattering by 
a r an do i p ■. im purity potential. For a single LL it was 
shownE3lljll3 that a delocalization transition can occur 
with increasing electron attraction and that repulsive in- 
teractions always lead to localization. On the other hand, 
for coupled LLs, a simple scaling analysis suggests that 
disorder would be irrelevant at least in the SLL phase.Q 
However, using a renormalization-group (RG) analysis, 
we show that disorder profoundly modifies the charac- 
teristic properties of these systems. It turns out that 
a delocalization transition persists in analogy to single 
LLs. Where Josephson inter-stripe couplings are irrele- 
vant, the delocalized phase can be identified with a new 
state of matter, which we dub disordered stripe metal 



(DSM). In contrast to the SLL state of the pure system, 
even in this delocalized phase, there exists only short 
ranged longitudinal CDW order due to impurity forward 
scattering. Because of this scattering process, we also 
find a strong tendency towards the destruction of trans- 
verse CDW order. Thus, the novel DSM state combines 
short ranged CDW order and quasi long-ranged longitu- 
dinal superconducting order with LL-like transport prop- 
erties. Interestingly, it has a much wider stability region 
in comparison to the pure system's SLL state. 



II. GENERAL MODEL AND 
RENORMALIZATION 



We assume a spin gaa in the LLs, as present in stripes 
in high-Tc compounds Jj The low-energy charge excita- 
tions of poninteracting stripes (labeled by j) can be 
describedE3 by the bosonic phase fields <I>j and their dual 
fields Qj with an action 

s° = ^Y. [ k(5xe,)2+t;^r(a,$,)2-2i9,$,a,e,- 

(1) 

The characteristic velocities vj^ and vj include forward 
scattering by intra-stripe interactions, whereas backward 
scattering is assumed to be irrelevaiii. 

Following Giamarchi and Schulz,t3 forward and back- 
ward scattering by weak impurities (denoted by IFS and 
IBS, respectively) is described in terms of the action 



^IBS 



V2 

TT 



/ Tij{x)d^^j, (2a) 
/ {G- (a;)e^^^*^ +li.c.}. (2b) 

A J XT 



r]j{x) and ^j{x) are Gaussian random variables with zero 
mean and correlations 77i(a;)?7j (x') = ^Djj5ijS{x — x') and 

£,i {x)£,j{x') = D^SijS{x — x'). a is the infinitesimal regu- 
larization length of bosonization. 
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In order to describe coupled arrays of stripes, we first 
include forward scattering due to density-density interac- 
tions between the stripes, such as screened Coulomb and 
electron-phonon interactions. The corresponding action 
reads 



5^ = ^y 



(3) 



In priiicjiiJe, analogous couplings between dxOj can be 
added,EHji!l but are dropped for simplicity. Their in- 
clusion can be achieved by an obvious generalization 
of our analysis, and will not modify our main results. 
Since is bilinear in $i, it can be included in the 
momentum-space representation of the action 5*" by re- 
placing vn in Eq. (|l|) by the q± dependent velocity 
VN{g±) = vn + Vq^ where Vq^^ is the Fourier transform 
of Vi (we use the spacing between stripes as transverse 
length unit). The restriction to low energy excitations 
underlying bosonization obviously requires the stability 
condition + Vq^^ > 0. The conventionally defined Lut- 
tinger parameter K — a/wj/uat is then generalized to 



Kiq±) = V^J/M<1±) = K/^l + Vq^lvN (4) 

in order to include inter-stripe forward scattering. 

Besides these forward scattering processes, we allow 
also for pairwise hopping between stripes, given by the 
transverse CDW and SC couplings 



S 



CDW 



cos 



2($. - $^)J , (5a) 
-&,)]. (5b) 



Since we assume the presence of a spin gap, HfWgle elec- 
tron hopping is irrelevant and can be ignoredO'EJ 

For the pure system without S^^^ and S^^^ , a spe- 
cific interaction Vi in Eq. (||) can render the CDW and 
SC couplings irrelevant in an intermediate region of K , 
leading to the SLL phase found in Refs. [^|j6|. The inter- 
action must be sufficiently strong and has to include at 
least nearest and next-nearest neighbors. 

In the presence of disorder, the scattering off impurities 
has to be taken into account. Let us first focus on the ef- 
fect of impurity forward scattering (|2a| ) in the absence of 
any inter-stripe couplings of type (H). This process then 
changes the SLL phase as described by 5*° -f into the 
disordered stripe metal (DSM). Introducing replicated 
fields and averaging over disorder still lead to a bilinear 
action with correlations 



uj'^ /vj + VN{qi_)ql vjj{qi_)q^^ 

n6''''vN{q±)/vj 
uj^/vj + VN{qi_)ql 



(6a) 



(6b) 



with q = {uj,q\\,qi_) and upper indices a, h as replica 
labels. 

We now examine the relevance of CDW and SC cou- 
plings and of impurity backward scattering (IBS) with 
respect to the DSM state using an RG analysis similar 
to that of Rcf. |l|Jl|. To first order in D^, C„i, and J„, 
we obtain the RG flow equations 



If 

"df 
dK 

~dr 

dCrn 

"df 



(3- 

Dr, 



2^3 



(7a) 
(7b) 
(7c) 

(7d) 
(7e) 



The scaling dimensions are 



with r — w/v is not renormalized to first order. 



a: 



.IBS 



CDW 



Kiq±), 



[1 — cos(r7igj^)] 



2a„ 



nAd%{qj_) 



A^^= / [l-cos(mgj.)]A--i(g^), 
with longitudinal momentum cutoff A ~ 1/a. 



(8b) 
(8c) 



III. PHASE DIAGRAMS 

Before we analyze specific models for the inter-stripe 
interaction, we discuss the general picture emerging from 
renormalization group. 



A. General Topology 

In the absence of disorder (13^ — ~ Q), K preserves 
its unrenormalized value y^vj/vN- Then the scaling di- 
mensions (|^) reproduce the expressions given in Ref . . 
For weak inter-stripe interactions V^^ vn the system 
is in the SC phase for X > 1, whereas it is in the CDW 
phase for K < 1. For a suitable choice of a strong inter- 
action Vq^ , an intermediate range of K can exist where 
all SC and CDW couplings are irrelevant and the system 
is in the SLL phase. The stability of the SLL phase is de- 
termined by the conditions K > if Cdw _ 2/ minm{c^} 
and K < K^'~^ = min„i{c~ }/2. Hereby we define 



c± = / [1 -cosTOgj_] ) 



±1/2 



(9) 



In the presence of disorder^ the strength of impurity 
forward scattering Z?^ increases exponentially under the 
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RG flow. This has two important consequences. First, 
the CDW order along the stripes becomes now short 
ranged as can be easily seen from the second term in 
Eq. (pal ). Second, it implies an exponential increase of 
^m^* for all m, i.e., the irrelevance of weak CDW cou- 
plings. Thus, impurity scattering transforms the SLL and 
CDW phases of the pure system into different phases. If 
impurity backward scattering (IBS) is irrelevant -i-this 
is the case in the entire stability region of the SLlJ3 - a 
novel phase is present which we call DSM phase. Un- 
like for the SLL, the stability of the DSM for smaU K 
is no longer limited by the CDW couplings but by IBS. 
Its phase boundary is determined by the relevance of SC 
couplings at large K and the relevance of IBS at small K. 
IBS leads to localization for a bare IBS strength larger 
than a critical value -D^,c (that depends on intra- and 
inter-stripe interactions). In this case 13^ diverges and 
K goes to zero under renormalization. For < D^ c, 
the system is delocalized, — * and K saturates at a 
finite value K* = K*{D^). For K below the critical value 
Kc = infinitesimal disorder produces localization 

(Dj_c = 0). For K > Kc, the system remains delocalized 
at finite disorder strength, < < D^ ^- 

To determine the phase boundary D^^^c of the localiza- 
tion transition for K > we integrate the flow Eqs. 
( [Za] ) and (|7c| ). In terms of the dimensionless disorder 
strength V — D^/{tt^Av]^) we obtain 



C C Kq 



with 



VN 
VN 



(10) 



(11) 



The critical disorder strength then follows from the con- 
dition that ^{K) = at its minimum a.t K — 3/c^ ^ 
Kc-. 

Vo.c = ^{Ko-Kc) + -ln^. (12) 
c c Ko 

In the delocalized phase, the renormalized value K* of 
K can be obtained from Eq. 1^ with V{K*) ^ 0. The 
phase boundary at large K between DSM and SC phase 



is given by the condition K* < K 



SC 



.{c™}/2, 



and the boundary at small K between the DSM and the 
localized phase is described by Eq. (12). 

The actual form of the phase diagram and, more im- 
portantly, the stability range of the SLL or DSM phase 
depends on the interaction Vq_^^ under consideration. To 
be specific, we will consider in the following two models 
for this interaction. 



B. Model A 

A minimal model that renders simultaneously all SC- 
and CDW-couplings irrelevant for some parameter region 
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FIG. 1: Phase diagram of model A for A = 10~^ in the 
absence of disorder. For large K, the system forms a 2D 
superconductor, while it is a 2D CDW for small k. The inter- 
mediate SLL exists only for suitable values of the parameter 
go. 



was suggested by Vishwanath et ala It assumes an inter- 
stripe interaction forward scattering leading to 



k{q±_) = k/[1 + Ai cos(gf_L) + A2 cos(2g_L)] 



(13) 



Within this model, the three parameters k, Ai and A2 
implicitly determine the intra- and inter-stripe interac- 
tion. The corresponding Luttingcr parameter is given by 
K"^ — K^/(l -I- (Aj -f A2)/2), and the inter-stripe potential 
Vi has a range of four stripe spacings. In analogy to Ref. 
^, we rewrite Ai and A2 as 



Ai 
A2 



4(1 - A) cos 90 
1 + 2 cos^ go 
1- A 



1 + 2 cos2 



90 



(14a) 
(14b) 



such that K[qi^ takes its maximal value k/A at q±^ = qq. 
By adjusting A to small values, CDW couplings are sup- 
pressed, while small k gives negative scaling dimensions 
for^SC-couplings, as pointed out by Mukhopadhyay et 
al.U One thus indeed finds regions in (go, k, A)-phase 
space, where the system is stable against all inter-stripe 
couplings of type ^ and is thus in the SLL phase. This 
is demonstrated in Fig. U which may be compared to 
similar plots in Refs. [plR . For sufficiently small A, win- 
dows of qo exist where the system evolves from a phase 
coherent SC through the 2D metallic SLL to a charge 
ordered CDW state with decreasing parameter k. Here 
and in the following, when both SC and the CDW cou- 
plings compete, the one that is most strongly relevant 
is assumed to determine the phase. The actual bound- 
ary between two such strong coupling phases might differ 
within a narrow corridor. Note however, that boundaries 
to either the DSM or SLL phase are obtained also quan- 
titatively correctly. 
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FIG. 2; Phase diagram for model A with A = 10~^ in the 
presence of disorder for two disorder values. The comparison 
of the phase boundaries at infinitesimal disorder T) = Q^ and 
finite disorder O = 5 • 10"^^ reflects the disorder induced 
renormalization of k. 



Now disorder is added while all other parameters are 
unchanged. The CDW and the SLL phases of the pure 
system become indistinguishable and merge to the metal- 
lic, short-range CDW-ordered 'disordered stripe metal' 
(DSM). Backscattering off impurities leads to localiza- 
tion in a large portion of the former CDW phase. The 
SC phase shrinks through downward renormalization of 
K by disorder (note that k, and K differ by a factor that 
is not renormalized) . In Fig. the boundaries between 
the three phases are given both for infinitesimal and for 
finite disorder. The latter shifts the boundaries to larger 

A cut through Figs. |l]and|^ at fixed go but with varying 
disorder is shown in Fig. |[ The SLL and CDW phases 
exist only for I? = 0. A delocalized phase can exist - due 
to inter-stripe forward scattering - even for purely repul- 
sive interactions (for example, go = O.SStt, A = 10~^ and 
K < 1.42 corresponds to repulsive on-stripe and repulsive 
inter-stripe interactions), as opposed to the strictly one 
dimensional electron gas with delocalization for i^T > 3 
corresponding to strongly attractive interactions. How- 
ever, the inter-stripe interactions corresponding to the 
values of A, k, go where the SLL or DSM exist may not 
be very realistic because of their strength. 



C. Model B 

Model A is constructed specifically in a way such that 
the inter-stripe forward scattering interactions (H) give 
rise to a nonmonotonous K[qi_) which allows for the si- 
multaneous irrelevance of CDW and SC couplings in the 
absence of disorder. For a large range of parameters go 
and A, this potential has oscillatory character in real 
space, which also may not be very realistic. 




0.035 0.1 0.2 



FIG. 3: Phase diagram in the plane spanned by disorder and 
interaction parameter k for fixed inter-stripe interactions with 
go ~ O.SStt and A = 10~^. With increasing disorder strength 
D, the boundaries of the DSM phase move to larger n keeping 
a nearly constant distance. 

A physically motivated choice for a potential that is 
monotonous both in real and Fourier space may be the 
screened Coulomb potential 

V{r) = -e-^"^ (15) 
r 

which we consider as model B. In Fourier space this 
model reads 

Vq^ = -Aln [1 + e-2M _ 2e"^cosg_L] . (16) 

Due to the stability condition Vq^/vM > —1, see Eq. 
(^), there is a critical amplitude Ac{fi), above which the 
model breaks down. 

Fig. |4| displays the stability of the model with respect 
to weak inter-stripe CDW and SC couplings in the ab- 
sence of disorder. No SLL is found, the system shows for 
all fj, and A a direct transition from the SC to the CDW 
phase for decreasing K. The addition of disorder leads to 
the phase diagram in Fig. |[ As opposed to model A, im- 
purity backscattering completely covers the CDW phase 
and thus leaves only two phases, the localized one and 
the SC phase. In contrast to the competition between 
CDW and SC couplings, the localization boundary is not 
given by the most relevant bare coupling. Since even 
weakly relevant IBS renormalizes K to small values, the 
disorder scaling dimension decreases while the Josephson 
coupling, provided a small enough bare value, ultimately 
becomes irrelevant, see Eqs. (0). Hence the boundary is 
given by the onset of relevance of IBS with respect to the 
pure system. It moves to larger K for increasing disor- 
der. No DSM phase is found now. Formally, the absence 
of a minimum of Vq^ inside the interval (0, tt) makes up 
for this latter qualitative difference in models A and B. 
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From a linear-response calculation we obtain the LL-like 
conductivity (which also determines the conduct ancella) 



2e2 



TTh (w + iO+)^/v* - VN{q±)q'^ 



(17) 
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FIG. 4: Phase diagram for model B for A — 0.78vm without 
disorder. 



0.8 



0.6 



0.4 



0.2 




model B 



localized 



4 




0.1 



0.5 



1.5 



K- 



FIG. 5: Phase diagram for model B for two disorder val- 
ues © = 0"*" and finite disorder D = 10~^ and infinitesimal 
Josephson coupling ^ = 0^; A = 0.78ujv; 



D. Correlations 

Having established the generic topology of phase dia- 
grams, we now address the nature of the possible phases. 
First we consider the delocalized DSM phase. It is de- 
scribed by the correlations (||) with the bare IFS am- 
plitude and the renormalized effective K* (in gen- 
eral, one has to use the renormalized but unrescaled 
quantities). We find a linear growth of the fluctua- 
tions of $ with longitudinal system size L, ($|(a;, t))^ — 

{c/2n'^)v]^'^DrfL, which leads to short-ranged longitudi- 
nal CDW correlations like for single LLs. On the other 
hand, IFS does not affect the quasi long-ranged longi- 



tudinal superconducting order (fluctuations of Q) of the 
pure system. Equally, the conductivity along the stripes 
is not affected by IFS since rij{x) is time independent .Ej 



representing a longitudinal metal. Here, ui represents a 
real frequency in contrast to Matsubara frequencies in 



Eqs. 



Notice that 



— vnK* . Since K* jumps 



from a finite value to zero at the localization transition, 
(T(q) behaves discontinuously there. In the transverse 
direction, CDW and SC correlations will be short-ranged 
since the corresponding couplings are irrelevant. In the 
presence of a spin gap (which suppresses single particle 
hopping) the irrelevance of the couplings also signals that 
the DSM is a transverse insulator. 

The localized phase is less amenable to an analytic de- 
scription since the divergence of IBS would necessitate 
a strong-coupling analysis. However, if the SC coupling 
is irrelevant, the localization transition and the localized 
phase share the qualitative properties of their ID coun- 
terparts. The inter-stripe couplings will lead to merely 
quantitative renormalization effects. Whether the local- 
ized phase is a random antiferromagnet or a pinned CDW 
depends on the mechanism generating the spin gap. The 
longitudinal localization length Lioc can be estimated an- 
alytically when the transition line is approached from the 
localized side. For K < K^, Ljoc ^ j^-i-p-^/i^i^-K/K,)) ^ 
For K > Kc, L\oc — A^^ exp{c/^/c(Vo^^V([c)) with c a 
numerical factor of order unity. Thus, the inter-stripe in- 
teractions influence the localizatiorLlength quantitatively, 
but the qualitative behavior foundEj for the localization 
transition in a single LL persists. 



IV. DISCUSSION AND SUMMARY 

In summary, we have examined impurity effects in ar- 
rays of coupled LLs. The competition between impurity 
backscattering, CDW and SC couplings allows for three 
different phases: a localized phase, a superconducting 
phase, and the disordered stripe metal. The latter two 
phases are delocalized since IBS scattering is irrelevant. 
While for a single stripe delocalization occurs only for 
strongly attractive on-stripe interactions {K > 3), for 
a coupled stripe array delocalization is possible also for 
purely repulsive on-stripe and inter-stripe interactions 
(forward scattering). 

The delocalized DSM phase is metallic in longitudi- 
nal direction and insulating in transverse directions. Its 
correlations for CDW order are short-ranged in all di- 
rections whereas superconducting correlations are quasi 
long-ranged along the stripes and short-ranged in the 
transversal direction. These experimentally accessible 
features should allow to identify the disordered stripe 
metal and to distinguish it from the SLL phase. 

In the above analysis we have determined the phase di- 
agram from a stability analysis of a Gaussian fixed point 
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- representing the stripe array with forward scattering by 
interactions and impurities ~ with respect to CDW and 
SC couphngs as well as impurity backward scattering. 
This stability analysis, reflected by the flow equations (Q) 
which are linear in Dj, Cm, and J7in, requires the weak- 
ness of these couplings. In principle, this approach does 
not cover strong-coupling phenomena: sufflciently strong 
couplings might drive transitions into SC or CDW phases 
which are less susceptible to disorder than the SLL. 

Although we cannot consistently access such strong 
coupling phenomena via our flow equations, they nev- 
ertheless can be used to determine crossovers related to 
the relative strength of non-Gaussian couplings. Since all 
CDW couplings are irrelevant at the DSM fixed point, 
we raise the question of whether a CDW phase can be 
reestablished if CDW couplings are sufficiently strong in 
comparison to disorder. We focus on the region where a 
CDW coupling Cm is relevant in the absence of disorder 
and where IBS is irrelevant (i.e., 5=2 — A^^^ > and 
^IBS ^ 3 ^j^g bare parameters). 

Irrespectively of the relative strength of CDW cou- 
plings and IBS, the presence of disorder implies a con- 
tinuous growth of Djj and thus also of A^'^^ which im- 
plies the irrelevance of Cm only on sufficiently large scales. 
Thus the question is, whether increases fast enough 
to achieve A^'^^ > 2 before a strong-CDW coupling 
regime is entered. This regime is entered when the di- 



mensionless coupling Cm = Ci/ {K^vn) becomes of order 
unity under renormalization before the disorder contri- 



bution to A 



CDW 



becomes of order 5. For weak Vn , or 



\, this is the case if 



Cm > 



2cDn 



(18) 



For CDWs (quasi-)long-ranged charge correlations can 
exist in the presence of disocder only in _D > 2 dimen- 
sions like for vortex latticesH3 Then the fermions would 
form a pinned (localized) Wigner crystal. However, in 
D — 2, the formation of CDW,order is prohibited by the 
proliferation of dislocations,Lll which ultimately render 
the CDW coupling irrelevant on sufficiently large scales. 
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